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The dynamics of transient patterns formed by front propagation in extended nonequilibrium 
systems is considered. Under certain circumstances, the state left behind a front propagating into 
an unstable homogeneous state can be an unstable periodic pattern. It is found by a numerical 
solution of a model of the Freedericksz transition in nematic liquid crystals that the mechanism 
of decay of such periodic unstable states is the propagation of a second front which replaces the 
unstable pattern by a another unstable periodic state with larger wavelength. The speed of this 
second front and the periodicity of the new state are analytically calculated with a generalization 
of the marginal stability formalism suited to the study of front propagation into periodic unstable 
states. 
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I. INTRODUCTION 

Front propagation is a common phenomenon in systems driven out from equilibrium [1]. It can occur when different 
regions of an extended system are prepared into unequivalent states. The interface between both regions is often sharp 
and moves at a characteristic velocity. In potential systems, for which the governing equations can be derived from a 
Lyapunov functional, the relative stability of the different states determines the direction in which the front moves. 
In nonpotential systems there are different stability criteria from which different kinds of fronts and front bifurcations 
can be obtained [2-4]. Interesting cases of patterns formed by chemical front propagation have been recently reported 

[5]- , _ ' 

Early rigorous work in front propagation addressed the advance of a uniform stable state into a uniform unstable 
one [6,7]. An important step forward was the recognition [8,9] that the rigorous results for the velocity of the front 
could be recovered from a simple marginal stability criterion [10-12]. It basically consists in studying the stability 
of the moving front neasecurityr its leading edge and stating that the selected velocity is the smallest among the 
ones for which the front is linearly stable. The physics justifying such prescription is that in the usual situation the 
state of the system accelerates along the manifold of unstable fronts until the first stable one is reached [13]. The 
marginal stability criterion cannot be applied when a nonlinear front dominates over the linear one [12,13], but it 
gives correct answers in many circumstances including nonpotential situations. The relation of the criterion with a 
structural principle has been recently addressed [14]. 

To apply the marginal stability prescription for front propagation into an homogeneous state, one writes the solution 
ip{x,i) of the relevant equations as ip[x,t) ~ exp(iqx + ui(q)t), linearize in t/>, and find the dispersion relation ui(q). 
The marginally stable front velocity v* is obtained from the simultaneous solution (q* , v*) of 
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(1.1) 



where q and q* should be considered as complex. This criterion has been successfully applied to situations where no 
exact results were available, such as vortex-front propagation in rotating Couette- Taylor flow [15,16] and Rayleigh- 
Benard convection [17]. Note that in these cases Re[</*] ^ and the stable state that advances into the unstable one 
is not homogeneous, but periodic. In this situation, front propagation has additional importance since it provides a 
natural mechanism for pattern selection. The marginal stability hypothesis, with some additional assumptions [8], 
gives the wavelength A of the pattern left behind the front: 
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Front propagation can also produce unstable periodic patterns. This was explicitly demonstrated in [18] for the 
extended Fisher-Kolmogorov model, an equation for which the only stable states are homogeneous. When this system 
is prepared in an initial condition with a sufficiently sharp interface separating the stable homogeneous state from an 
unstable homogeneous region, the advance of the front leads to the formation of a periodic pattern. This phenomenon 
could be considered as a generalization of the production of metastable phases by front propagation [19]. The situation 
is reminiscent of others in which a periodic pattern that is not stable appears after a linear instability [20]. At late 
times, the transient pattern evolves and decays. 

In the context of a model for the dynamics of nematic liquid crystals (which includes as a particular case the 
Cahn-Hilliard equation of spinodal decomposition) we find here a new process by which the unstable periodic pattern 
created by the propagating front approaches the homogeneous stable state: a second front advances through the 
periodic unstable pattern leaving behind a new periodic state with a larger wavelength. Under certain circumstances, 
the phenomenon can repeat itself until the homogenous state is attained. Our main goal in this paper is to study 
the characteristics of this second front and of the pattern created. We calculate the velocity of the second front and 
the periodicity of the pattern behind it. This is done by adapting the marginal stability theory to this new situation. 
The resulting equations are approximately solved and compared with numerical simulations. 

Front propagation into periodic unstable states has been discussed previously in [21] and [22]. In [21] the propagating 
front solutions of the 'amplitude equation' are considered with periodic unstable solutions as initial conditions. The 
difference with our situation is that in an amplitude equation description such initial condition can not be created by 
a primary front, since the front leaves a stable periodic pattern. In addition the decay of an initial unstable pattern 
is to a final stable pattern. More related to our case is the situation of [22], where complicated sequences of fronts 
and front splitting were found. A major technical difference between our situation and the ones discussed [21,22] 
is that the linear analysis around their stationary solutions lead to equations with constant coefficients, so that the 
usual formalism of marginal stability theory can be applied directly. We will deal here with linear equations with 
periodic coefficients, which is the generic case in the stability analysis around periodic solutions. Thus we need some 
adaptation of the marginal stability analysis to our case. 

The paper is organized as follows: In section II some general aspects of models leading to unstable patterns are 
reviewed. In section III, a marginal stability theory to characterize propagating fronts into an unstable periodic state 
is set-up. This theory is applied in section IV to a particular model realization for the Freedericksz transition in 
nematics. In section V, the velocity of the second front obtained from such theory and the wavelength of the pattern 
left behind are compared with values from the numerical solution of the model system, and the subsequent evolution 
of the pattern is discussed qualitatively from the observed simulations. Finally, section VI is devoted to conclusions. 



II. MODELS LEADING TO TRANSIENT PATTERNS 



In the general context of front propagation in potential systems, several cases can occur. The front separates states 
of different relative stability. The 'less stable' state can be either unstable or metastable and the 'more stable' can be 
again unstable or metastable or, obviously, stable. More important from the point of view of this paper, each state 
can be homogeneous or periodic. It is possible to find in the literature several examples of some of these cases. In 
the work of [18,11], unstable periodic patterns were shown to be formed behind a propagating front in the extended 
Fisher-Kolmogorov model (EFK). The EFK model is defined by the following equation for a scalar variable ip{x,t): 

1 j ; ( X; t) = d x 2 iP- 1 d x 4 iP + aP-bip 3 . (2.1) 

This is a generalization of the original Fisher-Kolmogorov (FK) equation (the model A of critical dynamics when a 
noise term is introduced) which is Eq.(2.1) with 7 = 0. 

In the context of the nematodynamic equations for liquid crystals, a model for the analysis of some general aspects 
of transient pattern evolution was analyzed in a recent paper [23]: 

i>(x,t) = (a-dZ)(dfy + c1>-b1> 3 ) (2.2) 

Eq. (2.2) defines a potential model which will be called Modified Cahn-Hilliard model (MCH). The original Cahn- 
Hilliard equation (CH) is the particular case of (2.2) with a = 0. The MCH model has also been considered by Salje 
[24] in a different context, for a partially conserved dynamics. It has also been shown to form unstable patterns by 
front propagation [25]. In the rest of the paper, we will assume a, 7, c, b > 0. When numerical values are needed, we 
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take c = 1 without loss of generality, and b = 3, as appropriate for modeling the Freedericksz transition in nematics 
with model MCH [23]. The free parameters of our study will be 7, and specially a. 

Let us summarize some of the properties of the four models defined so far (FK, EFK, MCH, and CH): All of them 
have in common the existence of at least three kinds of bounded stationary solutions: (a) tp(x) = ±^/c/& = ip±, 
(b)ip(x) = 0, and (c) a family of periodic solutions tpq(x) of fundamental wavenumber q in the range (0,q ma x) [23]. 
For the EFK model, q max = + 4:jc — l)/2j)~, and q max = yfc in the other models. For periodic boundary 

conditions this completes the list of stationary solutions for the MCH and FK models. There are additional solutions 
of the CH equation and EFK model not considered here. 

Regarding stability with respect to small perturbations, both solutions included in type (a) are linearly stable, and 
solution (b) is linearly unstable. All the periodic stationary solutions (c) are also unstable [23]. 

When a > c, there is a change of variables transforming the linear part of the MCH model into the linear part of 
the EFK. Since marginal stability theory makes use only of the linear portion of the evolution equations, the results 
obtained in [18] for the fronts in the EFK model can be immediately translated to the MCH model for a > c. The 
mapping is obtained by rescaling x and t in the MCH model into x' 
model is obtained in the new variables with 



•\Jaj{a - c) and t' = t/a. Then the EFK 



7 



(a-cy 



(2.3) 



Fronts advancing into the unstable uniform state in the EFK model leave behind one of the equilibrium homogeneous 
states ip± if 7 > 1/12, and a periodic pattern if j < 1/12 [18]. In the MCH model this last region is mapped into 
a < c(7 + 4v^) « 13.93c. 
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FIG. 1. Velocity of the first front calculated from (1.1) and wavenumber of the pattern left in the EFK model (a and b) and 
in the MCH model (d and e). Diamonds are simulation results for MCH model. 

For a < c, the EFK and MCH can not be mapped into each other and the marginal stability theory should be 
applied directly using the linear dispersion relation for the MCH model, oj(q) = (a + q 2 )(c — q 2 ), in Eq. (1.1). The 
procedure does not present any new difficulty. In Fig. 1 we show the velocity of the front moving into the homogeneous 
state and the wavenumber of the pattern left behind in the EFK model (as a function of 7) and in the MCH model 
(as a function of a with c = 1). For a < c, both the wavelength of the pattern and the velocity of the front in the 
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MCH model are small, hence a large amount of periods can be obtained with a small front velocity. This circumstance 
never happens in the EFK model. We also remark that it is precisely a < c the situation for which the MCH model 
gives results qualitatively different from the FK model [23] and appropriate to describe the Freedericksz transition in 
nematic liquid crystals. 

In Fig. 2 we show front propagation from a numerical solution of the MCH model for a = 0.2. Some algorithmic 
details follow: The integration of Eq.(2.2) has been performed by using a centered finite-difference scheme up to order 
(dx) 4 to approximate the spatial derivatives. To determine ip at t + dt a single step predictor-corrector method is used 
with dx = 0.25 and dt = 10 -4 [23]. The initial condition was ip(x,t = 0) = l/\/b for x < I and ip(x,t = 0) = for 
x > I, I being a suitable distance (for instance, / = 20). The boundary condition on the left was i]j{x < 0,t) = l/\/b. 
The size L of the system is initially large enough (L = 500) so that when the front arrives to the end, there is a well 
developed periodic structure on the left. Once the front arrives to the right end, the size of the system is slightly 
reduced by taking only the left region, and the remaining right part is taken as a new boundary condition, forced to 
be stationary from now on. By using larger systems, we have checked that the evolution at some distance from the 
right boundary is independent of this particular prescription for the boundary condition, as we expected from the 
local nature of the evolution discussed in [23]. 
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FIG. 2. Graph of xp(x, t) for t = 0, 20, 50, 125, 195, and 455 for a = 0.2. 

From the numerical solution we extract the following observations: First, the propagation of the front does not 
significantly affect the region x < I, so that this region remains in the ip + stable state. Second, a front advances into 
the initial homogenous state. Its velocity and the periodicity left behind are well determined by marginal stability 
theory (Fig. 1). Third and more important, a second front appears and moves into the periodic state left by the first 
front leaving behind a new periodic state. We interpret this new phenomenon in the following way: The advance of 
the first front creates in the system a periodic pattern that is very close to one of the periodic solutions of type (c) 
above. Since such solutions are unstable, they should decay. From the numerical results, the mechanism of decay is 
exactly the same as the decay of the initial ip = unstable state, that is, a front replaces the unstable state by some 
'more stable' state. This interpretation leads naturally to the study of stability properties and of marginally stable 
fronts propagating into periodic unstable states as relevant for the understanding of the second-front phenomenon. 

We have been unable to observe a second front in numerical simulations of the EFK model. We think however that 
this is a consequence of computational limitations and that the second front would appear at sufficiently long times 
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if the numerical scheme is accurate enough. This belief is supported by numerical simulations of EFK in which the 
initial condition is a periodic state of small wavelength (smaller than the one produced by propagation of the first 
front) connected to the homogeneous stable state. In this case a front was observed to propagate into the periodic 
state in very much the same way as in the MCH model. Thus we expect the general reasoning developed here to be 
of applicability in many transient-pattern forming systems. 



III. FRONTS PROPAGATING INTO PERIODIC UNSTABLE STATES 



There are several ways of formulating the marginal stability hypothesis. In all of them the dynamics of the front 
is analyzed in the leading edge, where the field is small enough so that the equation describing its evolution can be 
linearized. In the method of van Saarloos, the marginal speed theory is formulated in terms of the evolution equation 
of a profile with the form of a moving front [10-12]. A second approach is the pinch-point method developed in the 
context of plasma physics [26,27]. We will use here the method of the steepest descent (or saddle point) [8,13], since 
it is suitable to be generalized to front propagation into a periodic state, the aim of this section. 

We assume that the periodic pattern left by the first front is close to one of the stationary solutions with dominant 
wavenumber q = q%, that will be denoted by ip qi . We introduce ip(x,t) = ip qi {x) + 8ip(x,t) in the equation defining 
the model and linearize in Sip to obtain an expression of the form 

H= [C + U qr ]61> . (3.1) 

C is the linear operator giving the linear dispersion relationship u>(q) corresponding to linearization around the uniform 
solution ip = and U q stands for the remaining part: a periodic operator of periodicity q p related to qi. Eq. (3.1) is a 
linear equation with periodic coefficients whose formal solution is given by Bloch (or Floquet) theory [28,29] in terms 
of the eigenfunctions and eigenvalues of the linear operator C+U q . According to Bloch theorem, the eigenfunctions 
fk(x) of a periodic operator of period 2tt / q p are products of periodic functions of wavenumber q p and its harmonics, 
times plane waves of the wavenumber k which is used to label the eigenfunctions: 

oo 

fk(x) = e ik * *>r(k)e ir «'* . (3.2) 

r = — oo 

Here b r are complex amplitudes and A; is a continuous index in the infinite size limit. We assume that we are close 
to this limit, so that sums over h will be approximated by integrals. In addition to h, an additional discrete index, 
labeling eigenvalue branches, is required, but it will be omitted from our formulae to simplify the notation. The 
integral over h will always denote the sum over the values of h and over the additional discrete index. The Fourier 
transform of fk(x), denoted by fk(q) is given by: 

/oo 
dx e-'«* f k (x) = 2* ]T b r (k)S(q-k-rq p ) (3.3) 
r = — oo 

The integration is over the whole system. Given the initial perturbation Sipo{x) = Sip{x,t = 0), the solution of (3.1) 
reads in terms of these eigenfunctions and the corresponding eigenvalues e(k): 

/ + oo 

e c(k)t fk{x)C{k)dk , (3.4) 
- OO 

where we write the projection of the initial condition onto the Bloch basis as 

C(k) = j fl{x)84, {*)dx (3.5) 
fl is the complex conjugate of Using Parseval's theorem and (3.3), expression (3.5) can be simplified to: 

r. oo 

C(k) = (27T)" 1 / ft(q)SMq)dq = ]T 6j(A)^ (* + rq p ) . (3.6) 

r = — oo 

Using this result in Eq.(3.4) with Eq. (3.2) we obtain 
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oo P + OO 

SiP(x,t) = / dk br(k)bl l (k)6Mk + r'q p )e i(k+r ^>+< k ^ . (3.7) 

r,r' = -oo ■'- co 



Let us assume for a moment that only one of the terms in the sum (3.2) is relevant. This will be in fact the case 
for our particular model, but for completeness we will discuss the general case later. If only one term in the sums 
over r and r' is relevant, expression (3.7) is formally equivalent to the one used as the starting point in the steepest- 
descent presentation of marginal stability theory [8,13] for a front moving into an homogeneous state. The criteria for 
determining the marginal front velocity could then be generalized to our case with the front moving into the periodic 
state as follows. 

The interest is focused in the evolution of the perturbation Sip{x,t), which is supposed to evolve into a front, near 
its leading edge at large t. Restriction to the edge of the front validates the small Sip approximation performed in the 
linearization leading to (3.1). As the front moves in time, the study of Sip in the leading edge requires using a moving 
referential frame, that is, taking z = x — vt, where v is the velocity of that referential frame. The perturbation is then 
written as Sip(z,t) and this function has to be calculated at z ~ (leading edge) and t —> oo. Precisely for large t, 
the integral can be evaluated by the steepest descent method. 

With a single term r,r' = m taken into account and defining q' = k + mq p , Eq. (3.7) in the moving reference frame 
reduces to 

/+oo 
dq'\b m (q' -mq P )\ 2 6Mq')e h ^'^' z , (3.8) 
- oo 

we have also introduced h(q') = iq'v + e(k = q' — mq p ). 
According to the steepest descent method one has 

6i/>(z,t-> oo) ~ e h(q "^ t+iq " z (3.9) 

Subdominant time-dependencies and unimportant numerical constants have been omitted, q'* is a saddle point of 
the function h(q') analytically continued to the complex plane, that is dh(q')/dq' is zero at that point. Using the 
expression for h(q') one finds that this happens when 



de(k) 



dk 



= -iv (3.10) 

k = q'* —mq v 



This expression gives the position of the saddle point q'* in terms of the still undetermined velocity v of the reference 
frame. As it is assumed a priori that there exists a front moving at a constant speed, we state that when the reference 
frame moves at the same velocity of the front, v = v* , the perturbation Sip should remain finite and bounded in time 
in the vicinity of the leading edge (z ~ 0). From (3.9) this gives the additional requirement 

Re[%'*)] = (3.11) 

From this equation and (3.10) with v = v* one can determine the velocity of the front v* and the complex number 
q'* locating the saddle point. These equations are formally equivalent to (1.1), obtained for a front moving into an 
homogeneous state. The interpretation of q'* is the same as in such case: Eq.(3.9) shows that the real part of q'* 
gives the periodicity of Sip at the edge of the front, and its imaginary part characterize its steepness. The difference 
with the homogeneous case lies in the different eigenvalue spectrum e(k). 

Generalization of the above results to the case in which all the terms in the sum (3.2) are retained is now immediate. 
It is enough to estimate with the saddle point method each term in the sum (3.7). An equation identical to (3.10) is 
obtained for each value of r, which replaces m. There is now one value of q'* for each r, but (3.10) shows that all of 
them lead to the same k* . After some algebra we obtain 

6iP(x,t)^e< k "> t f k ,(x) , (3.12) 

where subdominant time dependencies and factors independent of x and t have been omitted. Using (3.2), and 
expressing (3.12) in the moving reference frame, z = x — vt, the condition of steady envelope of the front near the 
leading edge leads to the following set of equations (one for every r): 

Re[e(F) + i(k* + rq p )v] = (3.13) 

Since q p is a real number all these equations are in fact the same Eq.(3.11). Then (3.10) and (3.11) are valid also in 
this situation in which several modes in (3.2) are relevant. The only difference is that the state growing out of the 
unstable initial condition (see (3.12)) is not a single Fourier mode but a Bloch eigenfunction. 
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The wavelength A of the periodic pattern left behind by the moving front can be calculated following the prescription 
of Dee and Langer [8,10-12]. The main idea is to assume that the oscillations created at the leading edge by the 
linear instability will become quenched by the nonlinearities but their periodicity not modified. It is also necessary 
that one of the terms in Eq. (3.2) (say r = m) is the most important, so that Eq. (3.9) applies. In the moving frame 
of speed v* , linear theory predicts that the leading edge (3.9) oscillates at a frequency such that a number of nodes 
$ is created in the unit of time, with $ = 7r _1 (Im[e(fc*)] + i>*Re [</'*])• Behind the front, where the pattern has a 
wavelength A, the flux of nodes passing in the unit of time through a point fixed in the moving frame is 2v* /A. From 
this we get 

^(™ + ^'"l), (,.14, 

which is a generalization of (1.2) but now A is determined by the two different wavenumbers k* and q'* . The 
hypothesis behind this formula are that no nodes are created nor destroyed far from the leading edge, and that every 
node that linear analysis predicts to be created has to be really created. These assumptions work very well in the 
cases studied by van Saarloos and also in our case, but it should be remarked that they are not general and in some 
cases they seem to be violated, as for example in the case studied by Dee [21]. 



IV. SECONDARY FRONTS IN THE MCH MODEL 



The periodic pattern generated by the first front contains a dominant wavenumber qi. Our basic hypothesis to 
explain the second front phenomenon is that the pattern can be well represented by one of the stationary solutions 
of type (c) with such fundamental wavenumber. We study its stability with respect to new propagating fronts. To 
apply the formalism of the previous section, we first need the expression for ip q (x), the solution associated with a 
given wavenumber q. Then the spectrum of the operator C+U qp should be calculated and analytically prolongated to 
the complex plane. The solution of (3.10) and (3.11) gives then the velocity of the propagating front. None of these 
steps can be done exactly in an analytical way for the MCH model. They can in principle be fully done numerically, 
but since our aim is to explain some phenomena that have been already observed numerically, it would be preferable 
obtaining some approximate analytical expressions. We can find such approximation by working in the situation in 
which the wavenumber of the initial unstable configuration is close to yfc. In this case the amplitude of ip q is small, 
so that the associated U q can be considered a small perturbation to C. In addition, only the fundamental harmonic 
of q is important in this limit, a fact that simplifies calculations. By looking at Fig. 1, we see that the pattern created 
by the first front never has a periodicity qi = yfc (remember c = 1), but for sufficiently small a the approximation 
makes sense, since qi comes close to 1. 

To begin with, an approximate stationary solution of Eq.(2.2) (Duffing equation) is obtained by standard methods 
[30]: 

(4(c- g ?)/36) 1/2 S in( gi a:) (4.1) 

This solution is valid for qi < yfc, so that ip qi is small. Substitution of ip = ip qi + Sip into equation (2.2) gives 

= (a - d x 2 )(c + d x 2 - 3b rP 2 q .)6rP . (4.2) 

Comparing with (3.1), we identify U qp = — 3&(a — d x 2 )ip 2 t , which is small for qi < yfc and has a periodicity q p = 2qi. 
Writing 

^ 2 »= E A ' e2i,q ' x ( 4 - 3 ) 

(=-1,0,1 

qf)/(6b), and A\ = A_\ = — Ao/2, the eigenvalue problem for the equation (3.1) can be written as 

l 

e(k)bi(k) = uj(qi)bi(k) - r](qi) ^ A m bi_ m (k) ; / = -oo, • • • , +oo (4.4) 

m = — 1 

with oj(qi) = (a + qf){l — qf) and r/(qi) = 3b(a + q 2 ) , where qj = k + 2lqi. This expression represents, for every k, an 
infinite number of coupled equations for {&;}. 



where Aq = 2(c — 
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Several methods are known for calculating the spectrum of periodic operators approximately. Here we use the 
simplest weak coupling approximation in which U q is treated as a small perturbation of C [28]. The largest effect 
of the perturbation occurs at the boundary of the Brillouin zone k = ±qi where the degeneracy of the free spectrum 
uj(k) = <jj(k — q p /2) is removed creating a gap. However for k close to qi the two dominant free eigenfunctions are 
still eJ llqtX which correspond to r = , —1 in (3.2). Near the Brillouin boundary (4.4) can then be reduced to a 2 x 2 
matrix equation setting 6 r =0,r^0,— 1. The resulting equation have nontrivial solution for 

^(k) = Q^) (-A (u + «i) + (w + Wl ) 

± (w — wi) 2 + j4q(m 2 — uu\ + uf) + 2Ao(uiw — uw + uw\ — u\Wi) ^ , (4-5) 

where w = u>(k) , w\ = u>(k — 2g 8 ) , u = r)(k) and u\ = r)(k — 2g 8 ). It can be shown that e + (k) > for all k 
while e~ (k) is negative. So e + (k) gives the energy of the upper band an it is the required amplifying factor to be 
used in the integral (3.4). Expression (4.5) is only valid for k qi. In fact, if k = qi gives e + = 0, which is the exact 
eigenvalue associated to a uniform translation of the pattern. However far from k qi it lacks some properties of the 
exact spectrum, as for example the symmetry k — ► —k. 

The marginal speed of the front has to be obtained now by substitution of e + (k), with k taken to be complex, into 
Eq.(3.10) and (3.11). The resulting set of three real algebraic equations cannot be solved analytically. The method 
used to find the solution is as follows: the fourth order Taylor expansion of £ + (k) around a point k = k t is determined. 
Then, the derivatives involved in Eq.(3.10) are analytically calculated and the resulting set of equations numerically 
solved for v* and k* . As the Taylor expansion only fits the original function near a small neighborhood of k t , an 
iterative selfconsistent method is applied: k t is initially taken as the point where £ + (k) for k real takes its maximum 
value. A first approximation for v* and k* is numerically obtained and a new Taylor expansion oie + (k) is calculated 
around a new k t equal to the real part of the just obtained k* . The process is repeated until the solutions converge. 
The final solution turns out to be independent of the initial point k t provided that it is not too far from the maximum 
of e(k) with k real. Once the solutions v* and k* were obtained by this method, they were proved to verify the relevant 
set of equations. 

Within the approximation used to solve (4.4) it is found that as k moves away from the Brillouin boundary, &_i 
is the largest amplitude in (3.2) for the upper band of the spectrum (e + (k)). Then, as q'* = k* — 2qi is the most 
important wavenumber of the function fk*(x) the use of Eq. (3.14) gives the periodicity of the pattern left behind. 

In the next section, the velocity of the second front calculated from this theory as a function of qi and also 
the wavenumber of the pattern left behind will be presented and compared with numerical results from the direct 
integration of (2.2). 



V. NUMERICAL RESULTS AND COMPARISON WITH THEORY 



Figure 3 shows space-time plots of tl>(x,t) obtained from a numerical solution of (2.2) in gray levels and for several 
values of a. The horizontal axis represents space, and time runs along the vertical axis. White corresponds to regions 
with ip near the value of the stationary stable solution l/\/b and black to those with the value near -l/Vb. The 
location of the fronts is clear and the measured slopes determine their velocities. 

For the first front moving into the homogeneous phase, the velocity obtained fits the result obtained with the 
marginal stability theory (see Fig. 1). Each value of the parameter a determines through (3.14) a periodicity which 
is in agreement with the numerical solution. Such periodicity identifies a qi which is used as the initial wavenumber 
to apply the general results of IV for the second front. The values obtained from the theory and from the numerical 
solution for the second front are shown in Fig. 4a for several values of a. The mean periodicity left behind the second 
front is shown in Fig. 4b together with those obtained from Eq. (3.14). The error bars on the theoretical values are 
obtained when the value of qi assumed for the first front is modified by ±0.01 from the value predicted from marginal 
stability theory. This intends to take into account the fact that due to system size the wavenumbers in the numerical 
simulations are discrete and can then differ in this amount from the theoretical predictions. It is seen that this small 
indeterminacy on the value of qi can have quite a large effect in the final states. In addition this also suggest that 
small deviations from perfect periodicity in the state left by the first front can induce relevant changes in the state 
left by the second front. 
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FIG. 3. Space-time plots of ij)(x, t) in gray levels, for a = 0.0, 0.3, 0.6 and 0.9. Space is represented 
along the horizontal axis (system size: 137.5) and time runs along the vertical one from to 575. 
White and black correspond to values of xp(x,t) close to plus and minus \j\fh, the value of the 
homogeneous stationary solutions. 
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FIG. 4. (a) Speed of the second front for different values of the parameter a, obtained by using 
the marginal stability theory (O) and obtained from the numerical simulations (*). The vertical 
bar is a tolerance band (see text), (b) Comparison of the dominant wavenumber left behind the 
second front. 



Fig. 4a and Fig. 4b indicate that, as expected, our approximate analytical calculation is accurate for small values of 
a. In fact, the periodicity is closest to yfc for a = and the weak coupling calculation of the spectrum e(k) is only 
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meaningful for < y/c. The validity of the marginal stability hypothesis for the general case studied in section III of 
a front propagating into a periodic unstable state, can be checked applying the result of section IV for values of for 
which (4.5) is justified. These values do not generally correspond to a pattern created by a first front for a particular 
value of a. To this end, we have studied the appearance of a front for a = 0, and starting from the exact stationary 
periodic solution, instead of letting it be created by the first front. The stationary state was generated numerically 
taking tp(x, 0) from (4.1), and then use (2.2) until ip becomes steady with the growth of the harmonics modes. Then, 
a small part of this solution is replaced by ip + and we observe the evolution of the resulting configuration. For > 0.9 
the pattern is already very unstable and numerically it decays by roll annihilation before the appearance of the second 
front. On the other hand, for < 0.6 the velocity of the second front tends to be so small that the computer time 
needed to observe it becomes prohibitive. For intermediate we show in Fig. 5, the velocity of the front and the 
wavenumber of the periodic pattern left behind the front (determined from the average wavelength of the pattern). 
The agreement between the theoretical values and the numerical ones is quite good, taking into account that the 
calculations of e(k) is accurate for — ► y/c = 1. 

The wavelength of the pattern left by the second front becomes larger as the value of a is increased. Our numerical 
results indicate the possibility of the divergence of the wavelength for a finite value of a, in the same way than the 
wavelength of the first front diverges for a < c(7 + 4v^3). However the finite size of our simulations makes it difficult 
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FIG. 5. Velocity of the second front and wavenumber of the structure created when the front 
propagates into stationary states ip qt (x) with different dominant wavenumber (a = 0). The point 
represented by * corresponds to the value of q t obtained in the wake of the first front, for a = 0. 



The velocity of the second front is near one order of magnitude smaller than that of the first front. This means 
that the region occupied by the periodic pattern left by the first front will become larger as time goes on. Hence, it is 
expected that the decay in this region, far enough from the boundary condition, will occur finally through the local 
mechanism of roll destruction analyzed in [23], instead than through the second front mechanism. This is enhanced by 
the fact that the pattern left by the first front is not exactly periodic. In fact marginal stability theory is well founded 
only when the perturbation around the initial unstable state is localized in space. This assures that the dependence 
on the initial condition present but not explicitly written in (3.9) and (3.12) plays no role in the selection process. The 
difference between the pattern created by the first front and the closest periodic stationary solution is not localized 
and the mode of decay at very long times can be different from the front mechanism. Hence, as observed, the second 
front can only advance a finite (yet undetermined) distance before the pattern left by the first front decays by a bulk 
mechanism. 

This mode of spontaneous decay of the pattern, that is, unaffected by the boundary, seems to be also the general 
way of evolution of the structure left by the second front. In fact this pattern has a Fourier spectrum broader than the 
produced by the first front, and even the local wavenumber is slightly different from place to place. However, a very 
interesting case has been observed for a = 0.2. The pattern left by the second front decays towards the homogeneous 
state through the appearance of a third front. This is shown in the space-time plot of Fig. 6. The root for this 
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behavior seems to be the fact that the wavenumber of the second structure is fairly twice that of the first and this 
leads to a second structure with a very well-defined periodicity. For other values of a where there was an integer ratio 
between both wavenumbers, the second pattern has a wavenumber so small that we have not been able to observe 
numerically the decay of the structure in a reasonable computer time. In fact, when the wavenumber of the pattern 
tends to zero the main driving force for decay is wall interaction which leads to a very slow logarithmic evolution. 




FIG. 6. Space-time plot of ip(x,t) as in Fig. (3) for a is 0.2 (system size: 125.7 ; time from to 
1600. The decay of the unstable state = to = l/\/b through a sequence of three fronts is 
clearly observed. 

VI. REMARKS AND OPEN QUESTIONS 

The main conclusion of our work is the existence of a second front mechanism for the decay of unstable periodic 
pattern created by a first front in the MCH and CH models. This second front is well explained by marginal stability 
theory around a periodic state. As commented at the end of section II, we expect the second front phenomena to 
be a general feature of transient pattern forming models. The question of the observability of this phenomena in 
real experiments is related to the same question for the first front: due to thermal and environmental noise it is 
experimentally difficult to maintain a sample in a unstable state long enough to observe the advance of propagating 
fronts. As a consequence there are many experiments on front propagation into metastable states, but examples of 
front propagation into unstable states are scarce [15,17]. Nevertheless, the transient patterns observed during the 
Freedericsz transition in nematics [31] appear to be quite long lived, so that this system which is generally described 
by (2.2) [23,20] is an obvious candidate to look for the second front. 

Our theoretical approach has been intentionally restricted here to the situation in which some analytical or semi- 
analytical results could be obtained. A number of open questions which could be addressed by a more detailed and 
extensive numerical study include the possible divergence of the wavelength of the pattern behind the second front 
at a critical value of a, the interplay between front propagation and bulk decay and the conditions for observability 
of third and subsequent fronts. 

Finally it should be stressed that marginal stability theory in its present state is just an heuristic prescription. It 
would be desirable to establish rigorously its range of applicability. Such rigorous analysis has been only performed 
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so far for parabolic equations [7]. 
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